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PRELIMINARY  ASSESSMENT  OF  NEWLY  AVAILABLE  DIGITAL 
SEISMIC  DATA  FROM  KAZAKHSTAN 
Paul  G.  Richards  and  Won-Young  Kim 


INTRODUCTION 

This  report  concerns  a  unique  seismological  observatory  that  has  operated  since  the  1960’s. 
Located  in  Kazakhstan,  for  over  25  years  it  has  been  recording  and  archiving  up  to  50  channels  of 
digital  data  each  day. 

This  project  began  for  us  in  the  spring  of  1990  at  the  annual  workshop  of  the  Incorporated 
Research  Institutions  for  Seismology  (the  IRIS  Consortium).  Seismologists  from  the  Soviet 
Union  attended  and  told  us  that  they  had  permission  to  describe  a  geophysical  observatory  at 
Borovoye,  near  Kokchetav  in  Northern  Kazakhstan,  which  had  operated  digital  seismometers 
since  1966.  They  said  there  was  a  large  tape  archive  and  that  much  of  the  data  were  of  excellent 
quality  and  would  probably  be  made  available  on  request  -  even  the  data  for  nuclear  explosions. 

In  about  December,  1990,  one  of  us  made  a  data  request  for  signals  from  about  30  nuclear 
explosions,  and  from  several  earthquakes  near  nuclear  test  sites. 

At  the  next  IRIS  workshop,  held  in  March,  1991,  we  received  a  nine-track  tape  that  had  43 
examples  of  three-component  digital  data,  all  recorded  at  Borovoye,  Northern  Kazakhstan. 

In  this  report,  we  describe  some  of  the  data  we  were  given,  what  we  have  learned  about  the 
instrument  responses,  and  some  of  our  preliminary  analyses. 

STATION  CHARACTERISTICS 

The  Borovoye  seismic  station  (BRV),  with  coordinates  53°  03'  29"  N; 

70°  16'  58"  E  (see  Figure  1),  is  located  on  a  pluton  about  200  km  across  with  an  ill-defined  but 
deep  Moho  and  monolithic  granites.  The  station  owes  its  existence  to  its  ability  to  pick  up  signals 
very  well  from  the  Nevada  Test  Site.  The  data  at  Borovoye  of  .seismic  signals  from  Nevada 
explosions  have  P  waves  two  or  three  times  larger  that  would  be  expected  on  the  basis  of  globally 
averaged  magnitudes  for  Nevada  explosions.  The  status  of  "observatory",  with  continuous  analog 
recording,  was  ticcorded  to  Borovoye  in  1965.  The  early  deployment  of  seismic  instruments  was 
in  a  tunnel  driven  horizontally  into  granites  rising  to  the  south  of  a  lake  about  4  km  across,  but 
now  all  seismometers  are  in  a  15-meter  deep  vault  with  three  chambers. 

There  are  three  digital  systems.  The  first  BRV  digital  seismic  system,  known  as  KOD 
began  recording  in  1966  and  operated  continuously  from  1967  to  1973.  It  is  based  on  three 
component,  short-period  seismometers,  and  is  important  as  one  of  the  few  digital  systems  in  the 
late  1960's  -  early  1970's.  From  the  perspective  of  the  standard  convention  on  ground  motion 
polarity,  in  which  motions  in  the  up,  north  and  east  directions  are  taken  as  positive,  we  note  that  all 
KOD  signals  have  the  opposite  convention.  The  second  BRV  digital  system,  known  as  STR,  has 
operated  from  March,  1973  to  the  present  day.  It  consists  of  two  separate  parts;  STR-SS,  which 
is  intended  mainly  for  low-gain  recording,  and  STR-TSG,  which  includes  six  long-period  and 
seven  short-period  Kimos  seismometers,  most  recorded  at  two  gain  levels  for  a  total  of  20  data 
channels.  The  highest  sensitivity  is  100,000  counts/micron,  based  on  a  short-period  Kimos  with  a 
special  magnet,  and  a  low-noise  amplifier.  The  third  BRV  digital  system,  known  as  ASST  and 
installed  in  April,  1990,  uses  broad-band  forced-balance  seismometers  and  is  still  under  test. 

Figures  2  and  3  show  instrument  responses  for  .several  instruments  of  the  KOD.  STR- 
TSG,  and  STR-SS  systems.  Table  I  shows  parameters  of  instruments  in  each  of  the  four 
systems.  Figure  4  shows  calibration  pulses  for  the  TSG  .system,  together  with  the  synthetic 
response  we  have  obtained  at  Lament  in  temis  of  poles  and  zeros. 

EXAMPLES  OF  EXPLOSION  DATA 

In  Figure  5  is  shown  a  detailed  l(x;ation  map  for  .several  nuclear  explosions  at  the  Shagan 
River  test  site.  East  Kazakhstan.  In  Figure  6  is  shown  a  set  of  seven  seismograms  for  these 
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events,  all  recorded  at  Borovoye  on  the  same  instrument,  showing  excellent  data  for  P/i  and  the 
other  main  regional  phases  Sn,  and  Rg.  Notice  that  even  though  these  events  are  separated  by 
no  more  than  twenty  kms,  there  is  considerable  variability  in  amplitude  of  the  regional  phases 
recorded  at  BRV,  over  a  distance  of  about  7(X)  km.  Figure  7  shows  the  same  data,  with  an 
expanded  time  scale.  We  comment  below  on  the  Lg  signals,  which  we  have  compared  with 
NORSAR  and  WMQ  for  these  same  events.  Note  also  the  appearance  of  the  first  trace  in  Figures  o 
and  7,  which  turns  out  to  be  a  recording  of  a  double  explosion. 

For  the  Joint  Verification  Explosion  conducted  at  the  Shagan  River  test  site  on  1988 
September  14,  we  have  several  channels  including  the  TSG  short  period  and  long  period  channels 
(see  Figure  8)  and  a  three-component  broad-band  record  (see  Figure  9). 

As  representative  examples  of  the  Borovoye  data  for  explosions  at  Novaya  Zemlya,  Figure 
10  shows  a  three-component  short- period  record  for  the  event  of  1988  December  4.  It  may  be 
seen  that  the  short-period  Lg  is  quite  weak  at  this  distance,  almost  2400  km.  Figure  1 1  shows  a 
three-component  long-period  record  for  the  Novaya  Zemlya  explosion  of  1988  May  7. 

EXAMPLES  OF  EARTHQUAKE  DATA 

When  we  had  the  opportunity  in  1990  to  make  a  data  request,  a  puzzle  from  1969  came  to 
mind.  Figures  12  and  13  were  published  20  years  ago,  and  show  (respectively)  an  my  Ms  diagram 
with  an  event  marked  as  KA  appeiiring  in  the  part  of  the  diagram  normally  associated  with 
explosions,  and  the  data  for  this  event  reponed  at  the  Worldwide  Standttrd  Station  at  Kabul.  For 
its  body  wave  magnitude,  this  event  of  1969  May  1  has  a  very  low  surface  wave  magnitude.  But 
there  is  a  problem  with  measuring  tlie  surface  wave  magnitude  because  there  was  a  concurrent 
large  eanhquake  in  Tonga  that  generated  large  surface  waves  all  around  the  wrrrld.  Even  at  Kabul, 
so  much  nearer  the  Kazakhstan  event,  the  surface  waves  from  Kazakhstan  are  small  and  their 
magnitude  is  impossible  to  measure  accurately.  Figure  14  shows  the  three-component  Borovoye 
digital  data,  recorded  on  the  KOD  systetn  at  a  distance  of  about  10°  from  the  1969  event.  Figure 
15  shows  the  same  data  with  an  expanded  time  scale.  Figure  14  is  slightly  clipped  for  the  Lg 
phase,  but  otherwise  it  is  excellent  data.  For  example  it  is  good  enough  to  analyze  particle  motion 
in  detail,  and  just  after  the  P  arrival  we  find  a  quite  strong  signal  arriving  off  azimuth  by  several 
degrees. 

.Another  interesting  Kazakhstan  earthquake  occurred  on  1976  March  20,  and  in  Figure  16 
u  e  show  the  Borovoye  short-period  digital  data.  In  Figure  17  is  a  comparison  of  this  earthquake 
with  the  signal  from  the  Joint  Verification  Explosion  of  1988  at  Shagan  River.  We  have 
incctrporated  a  polarity  reversal  of  the  earthquake  to  make  the  comparison,  and  it  is  clear  that  the 
first  arrival  of  the  earthquake  is  a  tlilatation.  This  earthquake  was  a  thnist,  and  all  teleseismic 
arrivals  showed  a  compression.  Figure  18  shows  the  same  data  with  an  expanded  time  scale. 

EXAMPLES  OF  DOUBLE  EXPLOSIONS 

It  has  long  been  known  that  dtnible  explosions  arc  sometimes  carried  out  at  the  Kazakhstan 
test  sites  t Sykes  and  Ruggi,  1989).  Figure  19  shows  an  example  from  December,  1972,  reccided 
teleseismicaliy  at  NORS.AR.  The  larger  F  wave  is  slightly  dipped,  and  for  purposes  of  estimating 
an  Lg  magnitude  of  the  larger  event,  it  is  necessary  to  make  a  correction  for  the  contributitvi  from 
the  smaller  event  (see  Ringdal.  1989). 

From  Borovoye,  we  ret|ue.sted  two  examples  of  double  explosions.  This  data  was 
included  w  ith  the  tape  we  received  in  March,  1991.  In  Figure  20  we  display  the  Borovoye  short- 
period  data,  which  in  each  ca.se  is  for  an  explosion  at  Degelen  Mountain  followed  about  10  seconds 
later  by  a  larger  explosion  at  the  Shagan  River  test  site.  The  upper  example  of  Figure  20  is  the 
same  double  explosion  from  1972  that  was  illustrated  in  Figure  19.  The  low  er  example  in  Figure 
20  is  the  double  explosion  mentioned  above,  in  Figures  6  and  7,  occurring  in  1978. 

Since  this  is  digital  data,  we  can  look  just  at  the  two  Degelen  signals  on  an  expanded  scale, 
and  compare  them.  Such  a  compan>on  is  shown  in  Figure  21.  The  instruments  responses  are 
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slightly  different  in  1972  and  1978,  and  a  correction  has  been  made  to  obtain  data  with  the  same 
response.  It  is  cletir  from  Figure  21,  with  its  excellent  agreement  of  the  two  wave  forms  recorded 
about  six  years  apart,  that  these  two  events  were  located  very  close  to  each  other. 

ANALYSIS  OF  LG  PRECISION 

Next,  we  describe  the  measurements  of  the  Lg  signals  recorded  at  Borovoye,  and  the 
stability  of  this  measurement. 

Figure  22  shows  the  work  station  display  for  our  measurements  of  RMS  Lg.  We  choose  a 
window-length  and  plot  the  RMS  value  (more  correctly,  the  logarithm  of  the  root  mean  square 
averaged  over  the  window)  as  this  window  is  moved  down  the  seismogram.  As  shown  in  the 
figure,  the  RMS  value  is  indicated  along  with  the  corresponding  group  velocity.  We  have  found  in 
practice  that  our  results  on  the  stability  of  Lg  do  not  change,  whether  we  use  a  few  tens  of  seconds 
or  more  than  a  minute  for  the  Lg  window.  We  had  six  events  we  could  study  that  were  recorded 
on  the  same  instrument  and  that  were  not  clipped. 

Figure  23  plots  our  results  for  the  five  explosions  which  were  big  enough  for  a 
measurement  to  be  made  at  NORSAR,  almost  4000  km  away.  Figure  24  shows  the  comparison 
for  all  six  of  the  events  recorded  at  Borovoye,  compared  with  RMS  Lg  recorded  at  the  station 
WMQ  of  the  Chinese  Digital  Seismograph  Network.  In  both  cases  we  find  that  the  scatter  of  the 
data  about  a  best  fitting  line  is  quite  small,  having  a  one  sigma  value  of  about  0.05  magnitude  units 
on  three  degrees  (Figure  23)  or  four  degrees  (Figure  24)  of  freedom. 

CONCLUSIONS  AND  RECOMMENDATIONS 


We  now  have  37  examples  of  digital  data  from  the  Borovoye  Geophysical  Observator\'  for 
nuclear  explosions  conducted  at  Shagan  River,  Novaya  Zemlya,  and  Lop  Nor.  At  least  two  of 
these  examples  are  double  explosions.  We  also  have  data  on  about  10  earthquakes  (some  regional, 
some  teleseismic).  We  are  pursuing  the  question  of  calibration,  for  all  the  many  instruments  for 
which  we  have  so  far  received  data.  In  the  prtKess  of  acquiring  this  data  we  have  been  given  much 
information  on  what  is  clearly  the  most  sophisticated  seismic  station  ever  operated  in  the  Soviet 
Union. 


If  the  research  community  in  the  United  States  deems  these  data  valuable,  then  there  will  be 
some  consideration  of  an  effort  to  copy  as  much  as  possible  of  the  whole  archive  at  Borovoye  into 
a  modem  format  for  u.se  by  seismologists  in  Russia  and  Kazakhstan,  as  well  as  the  United  States. 

Such  an  effort  may  usefully  be  thought  of  as  two  separate  steps.  First  to  copy  all  the  wide 
tapes  (we  believe  there  are  about  7,000  of  them,  holding  a  total  of  about  70  gigabytes)  to  a  mexiern 
recording  medium.  This  would  not  be  a  a  major  project,  but  some  consideration  would  have  to  be 
given  to  special  staff  support  at  BRV  for  the  duration  of  the  copying  in  addition  to  hardware. 
Second,  to  work  with  BRV  personnel  on  the  design  of  a  suitable  new  format  for  the  archive,  to 
incorporate  all  of  the  available  information  on  instmment  responses  (now  recorded  only  on  log 
books  held  at  the  observatory)  and  to  implement  the  new  format  with  appropriate  quality  controls. 
This  second  step  would  be  a  major  project  requiring  a  plan  with  milestones  and  an  agreed 
schedule.  It  might  take  several  years  to  execute. 


We  are  not  advocating  either  the  first  steps  or  the  two  steps  together  at  the  present  time. 

But  if  the  data  we  have  so  far  acquired,  supplemented  by  the  data  we  have  been  told  will  be 
forthcoming,  is  found  to  be  of  high  quality,  then  we  expect  to  be  making  the  case  for  salvaging  and 
modernizing  what  clearly  is  a  unique  digital  archive  at  Borovoye.  " . . - 
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Table  1 .  Instrument  Response  Characteristics  at  Borovoye  Seismic  Station. 


System 

Seismometer 

Channel 

T,(a) 

(s) 

Ds(b) 

S  (c) 

(unit/pm) 

(Hz) 

Rate('^> 

(sample/s) 

Channel 

# 

KOD 

SKM-3 

High-gain'" 

3.5 

0.71 

30(K) 

2.0 

33.3 

1,3,4 

Low- gain  (Z)* 

3(X) 

2.0 

33.3 

- 

SKD 

- 

30.0 

0.71 

- 

- 

- 

- 

STsR-SS 

SKM-3 

High-gain 

2.0 

0.5 

2000 

1.8 

41.7 

1,7,8,9 

Low-gain 

200 

1.8 

10.4 

6,8,9 

SKD 

High-gain 

25.0 

0.71 

5.0 

0.14 

5.2 

2,3,4 

Low-gain 

0.5 

0.14 

5.2 

- 

STsR-TSG 

KS 

High-gain 

1.5 

0.71 

4500 

2.37 

38.5 

7,8,9 

KSM 

High-gain 

1.5 

0.5 

100000 

1.43 

38.5 
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(a)  Ts  =  Seismometer  natural  period  in  second,  (b)  Ds  =  Seismometer  damping  constant  critical 

damping  =  1/V2,  (c)  Sm  =  Nominal  sensitivity  in  count/micron,  (d)  fp  =  Normalization  frequency 
where  nominal  sensitivity  is  measured,  (e)  Rate  =  Sampling  rate  in  samples/second,  Acually  this  is 
base  channel  not  necessarily  a  high-gain,  *  Only  vertical  component. 
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Locatiiin  map.  Borovoye  is  about  500  km  north-west  of  the  Semipalatinsk  test 
site  in  Kazakhstan. 
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time  (sec) 

Figure  4:  Comparison  of  actual  and  synthetic  Figure  5:  Detailed  location  map  for  selected 

calibration  pulses  (based  ^n  poles  and  zeroes),  nuclear  explosions  at  the  Shagan 

for  instruments  in  the  STR-TSG  system.  River  test  site. 


Digital  Seismograms  (channel  #1)  at  BRV  from  KTS 
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S)  and  SP  (TSG-KSM)  channels,  JVE2 


Biirovoye  records  of  the  Joint  Verification  Experiment  (explosion  of 


Group  Velocity  (km/sec) 


Novaya  Zemlya  event,  mb=5.6  TSG-DS  channels 


Group  Velocity  (km/sec) 

Long-period  record  at  BRV  of  Novaya  Zemlya  explosion  ot  1988  May 


r,  0 


*.  f ) 


'1  f' 


Figur^J.-  An  ttib/Mg  diagram  from  Landers  (1972),  based 

on  Mg  values  computed  by  the  method  of  Marshall 
and  Basham  (1972). 


f  1 


3.34km/sec  2. 


Event  KA  5-1-69  at  KBL  6k 

The  three  components  of  the  surface  waves  geiieialcd  by  the  event  K  A  as 
recorded  at  Kabul.  The  energy  arriving  with  a  grotip  vcloeily  of  3  •  34  km  s'  '  may 
be  either  Love  or  liiglier  mode  Rayleigh  waves. 

Figure  13;  WWSSN  data  from  ttie  evi'nt  of  1969  May  1. 
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Event;  1969/5/  1  4:0:8./ 
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^igure  14,  but  with  expanded  time  seal 
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JVE2  (9/14/88)  &  Earthquake  (3/20/76) 
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As  Figure  17,  but  with  expanded  time  scale. 
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NOKSAK 


12/10/72  4;26;57.8  mb=5.6  (DM)  Double  explOSiOHS 
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•  Ihe  .sh(  r  t- fK' r  i  I'd  BRV  record.*’  of  two  double  explosions. 
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Events  on  8/29/78  2:36:58.0  &  12/10/72  4:26:57.8  ©Deqelen  mb=5.2  &  5.6 
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Figure  21:  Comparison  of  two  Degelen  explosions,  recorded  at  BRV .  C.on.pare  tlus  displ.iv  witli  the 

beginning  of  the  top  tr.aces  of  lip.iire  h  and  7,  to  appreciate  the  merits  of  digital  recordin;t 
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